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Abstract 

In this paper are presented the results of Monte Carlo simulations on the capa- 
bility of the proposed NEMO-km 3 telescope to detect TeV muon neutrinos from 
Galactic microquasars. For each known microquasar we compute the number of 
detectable events, together with the atmospheric neutrino and muon background 
events. We also discuss the detector sensitivity to neutrino fluxes expected from 
known microquasars, optimizing the event selection also to reject the background; 
the number of events surviving the event selection are given. The best candidates 
are the steady microquasars SS433 and GX339-4 for which we estimate a sensitivity 
of about 5 • 1CT 11 erg/cm 2 s; the predicted fluxes are expected to be well above this 
sensitivity. For bursting microquasars the most interesting candidates are Cygnus 
X-3, GRO J1655-40 and XTE J1118+480: their analyses are more complicated be- 
cause of the stochastic nature of the bursts. 
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1 Introduction 



The realization of a km 3 scale detector for astrophysical high energy neutrinos 
is one of the most important scientific goals of astroparticle physics. Due to 
their small interaction cross section, high energy neutrinos are expected to 
allow us to observe inside dense environments in our Universe and to extend 
our knowledge beyond the distances that can be explored with 7-rays and 
cosmic ray observatories [1,2]. 

A number of astrophysical sources of high energy neutrinos, both galactic 
(SuperNova remnants, microquasars, ...) and extragalactic (Active Galactic 
Nuclei, Gamma Ray Bursts, ...), have been proposed. Since neutrino produc- 
tion requires a hadronic component inside the jet, neutrino detection will be 
a fundamental probe to investigate jet composition and particle acceleration 
processes taking place in these sources. 

The main aim of this paper is to study the sensitivity of the proposed NEMO- 
km 3 telescope [3] to high energy neutrinos emitted by known microquasars. 
In particular we simulated the interaction of microquasar muon-neutrinos in- 
side and in the proximity of the NEMO-km 3 neutrino telescope in order to 
estimate the expected muon event rate. In our calculations we accounted for 
the atmospheric neutrino and muon fluxes reaching the detector in order to 
optimize the event selection and reconstruction procedures. 



2 Neutrino fluxes from microquasars 

Microquasars are Galactic X-ray binary systems which exhibit relativistic jets, 
observed in the radio band. The presence of an accretion disc and relativistic 
jets make microquasars similar to small quasars (see ref. [4]). Several authors 
propose microquasar jets as sites of acceleration of charged particles up to 
energies of about 10 16 eV, and of high energy neutrino production. 

Levinson & Waxman [5] proposed a theoretical model for neutrino emission 
from microquasar jets, considering the production of neutrinos from the inter- 
actions of protons, accelerated in the jet, with synchrotron photons, emitted 
by electrons, or with external X-ray photons of the accretion disk. Accord- 
ing to this model, microquasars may be intense sources of high energy muon 
neutrinos, in the energy range between 1 and 100 TeV. Applying the Levin- 
son & Waxman model, Distefano et al. [6] estimated the expected neutrino 
fluxes from an ensemble of known microquasars, showing that several of these 
sources may be detected with the future km 3 neutrino telescopes. 
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Bednarek [7] considers the case of microquasars in which a Wolf-Rayet star 
supplies matter onto the compact object. Together with the scenario proposed 
by Levinson and Waxman, Bednarek considers also two regions where high en- 
ergy neutrinos may be produced. In the jet region II (see ref.[7]) neutrons re- 
sulting from photo-disintegration of accelerated nuclei produce neutrinos after 
interactions with the inner disk matter. In the jet region III the process takes 
place in the proximity of the Wolf-Rayet star surface. In particular Bednarek 
discusses the case of Cygnus X-3 and calculates the neutrino spectra (see Fig. 
2 in [7]) produced by a power-law spectrum of nuclei 7 A K , exploring different 
possible values for the spectral index k and for the Lorentz factor cut-off *y™ ax . 
Another component of the neutrino flux may arise from neutrons injected by 
mono-energetic nuclei having a Lorentz factor ^y™ n Rj 

10 5 . 



Aharonian et al. [8] discuss different possible scenarios for the production of 
the 7-ray flux recently observed from LS 5039 by the HESS Collaboration [9]. 
They consider both leptonic and hadronic production mechanisms and argue 
in favor of a TeV photon flux originating from pp interactions. They suggest 
therefore that the 7-ray flux should be accompanied by a TeV neutrino flux 
of about 1CT 12 cm -2 s -1 , or even a factor of 100 larger. 



In this paper we discuss the results obtained simulating neutrino fluxes calcu- 
lated with the three models reported above. 



Christiansen et al. [10] recently published predictions on the emission of high 
energy neutrinos from "windy microquasars" . Their calculations, based on the 
model proposed by Romero et al. [11], describe the emission of high energy 
neutrinos and 7-rays from pions created in the inelastic collisions between 
relativistic protons ejected by the compact object and ions in the stellar wind. 
As an example they consider microquasar LS I +61 303, estimating an average 
event rate of 3-5 TeV muon neutrinos per kilometer-square per year. This 
source is outside the field of view of a Mediterranean neutrino telescope and 
it is not considered in this work. 



Romero & Orellana [12] discuss the case of microquasars that show strongly 
misaligned jets with respect to the perpendicular to the orbital plane. If the 
donor star is an early-type star, the jet could collide with the stellar wind, 
producing 7-rays and neutrinos. The expected neutrino fluxes are in general 
too low to be detected by the planned km 3 telescopes and are not considered 
in this paper. 
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3 The NEMO project 



The NEMO Collaboration [13] is performing R&D towards the design and 
construction of the Mediterranean km 3 neutrino detector. The activity was 
mainly focused on the search and characterization of an optimal site for the 
detector installation and on the development of a feasibility study for the 
detector. 

After eight years of activities in seeking and monitoring marine sites in the 
central Mediterranean Sea, the Collaboration selected a large marine area 
(centered at the coordinates Lat. 36° 25' N, Long. 16° 00' E) located about 80 
km from the Southern cape of Sicily, Capo Passero, in the Ionian Sea plateau 
as an optimal site for the deployment of the km 3 telescope. The measured 
average values of the absorption length for blue light (440 nm) is ~ 70 m, 
close to the one of optically pure sea salt water; measurements also show that 
seasonal variations of the light absorption length are negligible [15]. The opti- 
cal background, measured at ~ 3000 m depth and for a 10" PMT at 0.5 s.p.e. 
threshold, has an average rate of 30 kHz, compatible with the noise expected 
from 40 K decay. An exhaustive report on the Capo Passero site properties is 
available in [15,16]. 

NEMO proposes a preliminary project for the km 3 detector based on present 
technological and budget constraints. The detector architecture consists of a 
square array of structures called towers, hosting the optical modules and the 
instrumentation. The tower is a three dimensional flexible structure composed 
of a sequence of "storeys" (which host the instrumentation) interlinked by a 
system of cables and anchored on the seabed. The structure is kept vertical by 
an appropriate buoyancy on the top. In its working position each storey will be 
rotated by 90°, with respect to the upper and lower adjacent ones, around the 
vertical axis of the tower. The final features of the tower (number and length 
of storeys, number of optical modules per storey, distance between the storeys, 
distance between the towers) is under study with the goal of optimizing the 
detector performance. 

The km 3 telescope, simulated in this paper, is a square array of 9 x 9 towers 
with a distance between towers of 140 m. In this configuration each tower 
hosts 72 PMTs (with a diameter of 10"), namely 5832 PMTs for the whole 
detector with a total geometrical volume of ~ 0.9 km 3 . We considered an 
18 storey tower; each storey is made of a 20 m long beam structure hosting 
two optical modules (one downlooking and one looking horizontally) at each 
end (four OMs per storey). The vertical distance between storeys is 40 m. A 
spacing of 150 m is added at the base of the tower, between the anchor and 
the lowermost storey. 
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As an intermediate step to ensure an adequate process of validation, NEMO 
built a demonstrator which includes most of the critical elements of the pro- 
posed km 3 detector: NEMO Phase-1 ([13] and ref. therein). It was realizated 
at the Underwater Test Site of the Laboratori Nazionali del Sud in Catania, 
where a 28 km electro optical cable, reaching the depth of 2000 m, allows the 
connection of deep sea intrumentation to a shore station. The NEMO Phase- 
1 system is composed of a junction box and a mini-tower (with 4 storeys). 
This allows to test the proposed mechanical solutions as well as the data 
transmission, the power distribution, the timing calibration and the acoustic 
positioning systems. NEMO Phase-1 is operating since December 2006 [14]. 

Although the Phase-1 project will provide a fundamental test of the technolo- 
gies proposed for the realization and installation of the detector, these must 
be finally validated at the depths needed for the km 3 detector. For these mo- 
tivations the realization of an infrastructure on the site of Capo Passero has 
been undertaken [13]. 

A further R&D program will also be developed by the Collaboration within 
the KM3NeT Design Study [17]. 



4 The simulation codes 



In this section we describe the Monte Carlo (MC) simulation tools and pro- 
cedures to study the detector response to high energy neutrinos. As a first 
step we generated samples of muon events having vertexes within an appro- 
priate volume of water. Three sets of muon samples were generated: atmo- 
spheric muons, muons induced by atmospheric neutrinos and muons induced 
by neutrinos coming from point-like astrophysical sources. These events were 
then propagated inside the detector, using the simulation tools developed by 
the ANTARES Collaboration [18] described in ref. [19]. The codes simulate 
the emission and the propagation of Cerenkov light radiated by muons and 
their secondary products, then record photo-electrons signals on PMTs. Subse- 
quently, background signals due to 40 K decay are generated in an appropriate 
time window and added to the event. Another code [19] is then used to recon- 
struct muon tracks. These codes, developed for detector configurations with 
smaller size and different geometries, were modified for a larger geometry, as 
described in [20,21]. A brief description of the simulation steps is given in the 
following. 
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4-1 Generation of neutrino-induced muon events at the detector 



Neutrino-induced muon events are generated using the event generation code 
developed by the ANTARES Collaboration [19]. The code generates neutrinos 
that interact inside and in the proximity of the detector, producing detectable 
muons. Interacting neutrinos are generated with a power law spectrum e~ x , 
where X is a generic spectral index. These events are weighted, as described 
below, in order to reproduce the expected spectrum. The direction of incident 
neutrinos coming from a point-like source is generated taking into account the 
change of the relative position of the astrophysical sources with respect to the 
detector rest frame, due to the Earth rotation. In the case of a diffuse neutrino 
flux, such as atmospheric neutrinos, the neutrino direction is generated with 
an isotropic distribution. 

Once the interacting neutrino is generated (i.e. neutrino energy e u , direction 
(i?^, tp v ), interaction vertex (x, y, z) and event weight W even t are assigned), the 
code simulates the neutrino Charged Current weak interaction, producing the 
muon kinematics. Since the detection efficiency of the detector is poor for 
muon energies lower than 100 GeV [21], the lowest neutrino energy considered 
in the simulations is 100 GeV. At these energies, the main contribution to the 
neutrino interaction, is from Deep Inelastic Scattering (DIS), here simulated 
using the code LEPTO [22]. Once the neutrino-induced muons are generated, 
they are propagated up to the surface of the detector sensitive volume. This 
volume is defined as a cylinder with size large enough that the Cerenkov light 
emitted outside its volume has a negligible probability to reach any PMT. 
In particular, the cylinder size extends beyond the instrumented volume by 
3 times the blue light absorption length (~ 70 m). The muon propagation is 
performed using the code MUSIC [23]. 

In the generation procedure a weight W event is assigned to each event in or- 
der to normalize the number of generated events to the expected neutrino 
flux. The event weight is calculated multiplying the expected neutrino spec- 
trum (dn v 1 1 de u dSdt) expected (evaluated at the generated neutrino energy and 
direction) times the generation weight W gen : 



W gen is defined as the inverse of the simulated neutrino spectrum. In the case 
of a diffuse neutrino flux, it is given by: 




expected 
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where X is the generation spectral index, Ie is the integral of the generation 
spectrum shape e~ x over the whole simulated neutrino energy range, 1$ is the 
integral of the solid angle in which events are generated. N total is the num- 
ber of simulated events and t gen is the event generation time. The neutrino 
interaction vertices are randomly generated within the volume V gen that com- 
pletely contains the detector sensitive volume and with size large enough that 
neutrinos interacting outside cannot produce detectable muons. In particular, 
its size extends beyond the detector sensitive volume by the maximum range 
covered by the neutrino- induced muons. {o'cc{ e v)pN a) -1 is the neutrino CC 
interaction length in a medium of density p, ccc{ e v) is the CC neutrino in- 
teraction cross section and Na is the Avogadro number. Eq.(2) takes also into 
account the neutrino absorption in the Earth, through the transmission prob- 
ability PEarth(£v,$v), which depends on neutrino energy and direction [24]. 
In the case of a point-like source neutrinos, the factor 1$ does not appear in 
relation 2. 



4-2 Generation of atmospheric muons 



Atmospheric muons are generated according to the Okada formula [25], which 
parameterizes the integral energy spectrum of muons, reaching a given depth 
D in seawater, as a function of the muon energy and direction. Since this 
parameterization is a complicated function of many parameters, a weighted 
generation procedure, similar to the one used for neutrino events, is used. 
Muon energies are generated with a given spectral index X and muon direc- 
tions are isotropically generated, in 27r solid angle (only downward directed). 
For each generated event a weight W event is calculated multiplying the Okada 
atmospheric muon spectrum times the generation weight W gen 



£ u lEl&t gen A g 

Ntotal 



Wgen — T-j , \0) 



where A geom is the detection geometrical e. the detector sensitive sur- 

face, projected on a plane perpendicular to the incident muon direction. The 
other parameters have the same meaning as in Eq.(2). In order to take into 
account the detector height (~ 1 km), the Okada spectrum is evaluated at the 
depth where the muon track intersects the detector sensitive volume. 



4-3 Muon propagation inside the detector 



The subsequent steps of the simulation are: the propagation of muons inside 
the detector, the generation of the Cerenkov photons and the simulation of the 
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PMT response. The muon propagation is performed with MUSIC [23] in order 
to take into account high energy radiative processes like bremsstrahlung. The 
Cerenkov light produced by secondary electrons is also taken into account. In 
the code we consider the light absorption length as a function of wavelength, 
according to the water optical properties measured in the Capo Passero site. 
Once the PMT hits are generated, spurious PMT hits, due to the under- 
water optical noise ( 40 K decays), are introduced. It is assumed that optical 
background produces uncorrelated single photo-electron (s.p.e.) signals in the 
PMTs. The PMT response function is taken into account. In the simulations 
made for this work the optical background signals are generated assuming 
an average rate of 30 kHz for 10" PMTs, corresponding to the average value 
measured in Capo Passero. 

4-4 Event reconstruction 

It is possible to reconstruct the direction of the muon track, using information 
from arrival times of photons in the PMTs, hit amplitudes (proportional to 
the number of photo-electrons in the PMTs) and PMT positions. The used 
reconstruction strategy aims at the rejection of hits due to background and 
selects a sub-set of hits compatible with the muon track. The track recon- 
struction algorithm used in this work is a robust track fitting procedure based 
on a maximum likelihood method [26]. 

In order to reduce the number of hits due to background the first step consists 
in the rejection of hits with amplitude smaller than 0.5_p.e. then a causality 
filter with respect to the highest amplitude hit is applied_J. The causality filter 
used in this work is based on the following condition: 

(\dt\ - dr/v < 20 ns) AND (||dt| - dr/c\ < 500 ns), (4) 

where \dt\ is the absolute value of the time delay between two hits, dr is 
the distance between the PMTs where the hits are detected, v is the group 
velocity of light in water. The first condition is based on the fact that the front 
propagation of direct Cerenkov photons moves in water with velocity v ~ c/n 
(n = 1.35). The second condition takes into account that, in the case of large 
distances, light absorption does not allow Cerenkov photons to propagate far 
from the initial muon track. The value of 500 ns takes into account the late 
Cerenkov photons, due to light scattering, 5-rays and showers. 

The reconstruction algorithm starts with a linear prefit applied on a sub-set of 

3 This is an important step since the number of noise hits may be comparable or 
larger than the number of muon hits in case of high noise rate or low muon energy. 
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hits that passed the causality criterion and in which there are at least three hits 
that satisfy the following trigger condition: the hit has an amplitude higher 
than 3.5 p.e. or the hit belongs to a "local coincidence". A local coincidence is 
defined as two or more simultaneous hits (At < 20 ns) in the two PMTs at the 
bar edge. Starting from the result of the pre- fit, a sequence of fit procedures 
using all the hits that passed the causality criterion is applied. A further hit 
rejection based on the track parameters is also applied. 



4-5 Criteria for atmospheric background rejection 



At the reconstruction level the fraction of background events is very high. In 
order to reduce it, an event selection is applied. The used reconstruction algo- 
rithm is a robust track fitting procedure based on a maximization likelihood 
method [26]. In this work, we used, as a goodness of fit criterion, the variable: 

A=-l^ + 0.1(7V comp -l), (5) 



where \og(C)/N DOF is the log-likelihood per degree of freedom (N DOF ) and 
N cornp is the number of solutions found by the reconstruction program (see ref. 
[26] for details). This quality cut is applied together with the selection criteria 
listed here: 

• the muon must be reconstructed with a Zenith angle ^ ec < $™ ax , in order 
to reject down-going events; 

• the variable A must be greater than a given value A cut ; 

• only events reconstructed in a circular sky region centered on the source 
position and having a radius of run are considered. 



5 NEMO-km 3 sensitivity to microquasars neutrino fluxes 



The detector sensitivity was calculated according to the Feldman and Cousins 
approach [29]. The 90% c.l. sensitivity to a neutrino flux coming from a point- 
like source is given by 

f _ ^90 W f th (a \ 

V 



where /x 90 (&) is the 90% c.l. average upper limit for an expected background 
(atmospheric neutrinos + muons) with known mean value b and no true signal 
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[29] , fl h is the theoretical neutrino flux from a given microquasar that induces 
a mean signal N™. During the calculation of the sensitivity, an event selection 
is applied as described in the previous section. The optimal values of $™ ax y 
A cut and r bin are chosen to optimize the detector sensitivity. 



5.1 Atmospheric neutrino background simulation 



In order to evaluate the contribution of the atmospheric neutrino background 
with respect to the signal produced by microquasars, we simulated a total 
number N to tai = 5.1 ■ 10 10 interacting atmospheric neutrinos, with a genera- 
tion spectral index X = 2. Neutrinos were generated in the energy range 100 
GeV - 10 s GeV with a Air isotropic angular distribution. The corresponding 
total number of reconstructed neutrino-induced muons is ~ 6.6 • 10 5 . The at- 
mospheric neutrino events are weighted to the sum of the Bartol flux [27] and 
of the rqpm model [28] to take into accunt the contribution of prompt neutri- 
nos. Considering an observation time of 1 year, the total weighted number of 
reconstructed atmospheric neutrinos (up-going + down-going) is ~ 4.7 ■ 10 4 
(quality cuts are not yet applied). 



5.2 Atmospheric muon background simulation 



Atmospheric muon events are also generated applying a weighted generation 
technique. We generated N tota i = 2.5 • 10 7 muons, in the energy range 1 TeV 
- 1 PeV, with a spectral index X = 3 and N total = 4 • 10 7 muons in the range 
100 GeV - 1 TeV, with a spectral index X — 1. Muons are generated with an 
isotropic distribution as described in sec. 4.2. The corresponding number of 
reconstructed events is ~ 9 • 10 5 in the range 100 GeV - 1 TeV and 3.1 • 10 6 in 
the range 1 TeV - 1 PeV (in total ~ 4 • 10 6 ). 

According to the Okada parameterization, the expected number of recon- 
structed muons is ~ 1.3 • 10 9 per year. The statistics of the generated muon 
sample corresponds to only a few days of data taking. Since reconstructed 
events have a flat distribution in Right Ascension (RA), it is possible to project 
the full sample of simulated events in a few degrees bin ARA. This allows to 
obtain statistics of atmospheric muon background corresponding to a time of 
~ 1 year of data taking for a region of a few degrees around each single source. 
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5.3 Simulation of neutrino events from microquasars 



For each source, a number N tota i = 10 9 of interacting neutrinos was simulated 
in the energy range 1-100 TeV in order to reproduce the neutrino fluxes 
calculated by Distefano et al. [6]. The generation spectral index is chosen 
to be X = 1, in order to guarantee a good event statistics at the highest 
energies. For microquasar LS 5039, N totat = 10 8 events were generated in the 
energy range 0.1-1 TeV to simulate neutrinos events according to the model 
proposed by Aharonian et al. [8]. For microquasar Cygnus X-3, N to tai = 10 8 
events were generated in the energy range 0.1-1 TeV and N to tai = 10 9 events 
with energies 100 - 1000 TeV to simulate Bednarek's theoretical fluxes [7]. 

Using the information on reconstructed muon directions we computed the 
equatorial coordinates (declination and right ascension) of the detected muons, 
and we mapped the reconstructed events in equatorial coordinates. This al- 
lows to obtain an image of the source and to evaluate the spread of the recon- 
structed muon directions (due to neutrino interactions and to the muon track 
reconstruction error) around the source coordinates. 

The 2D angular distribution of the reconstructed muons for SS433 and GX339- 
4 (the two persistent microquasars expected to produce the largest number of 
events) is plotted in Fig. 1. The angular resolution is calculated as the median 
of the angular spread distribution of the reconstructed tracks with respect to 
the source direction. It is 1.1° for SS433 and 0.9° for GX339-4 (see Fig. 1). 
The angular resolution value depends on the source theoretical spectral index 

(r = 2). 



SS433 GX339-4 




289 288 287 257 256 255 254 

RA (deg) RA (deg) 



Fig. 1. Weighted MC event distributions for the microquasars SS433 and GX339-4, 
referring to 1 year of data taking. Both maps are obtained from 2D histograms with 
0.06° x 0.06° bin size. The circle has a radius equal to the median of the angular 
spread distribution of the reconstructed tracks. It represents the region where 50% 
of the events are reconstructed. 



12 



5.4 Detector sensitivity to neutrino microquasars 



The detector sensitivity was calculated for a livetime of 1 year, simulating a 
neutrino flux with spectral index T = 2 in the energy range 1 - 100 TeV. The 
study was carried out for each microquasar, since the sensitivity is a function 
of the source astronomical declination. 

Tab. 1 gives the detector sensitivity /^go for each microquasar and the corre- 
sponding values of $™ ax ; A cut and r^ n which optimize f u ,9o- Fig. 2 shows the 
detector sensitivity for the studied microquasars as a function of the decli- 
nation; the sensitivity flux limit increases with increasing declination, due to 
the decrease of the time per day spent by the source below the Astronomical 
Horizon (with respect to the latitude of the Capo Passero site). In Tab. 1, 
we present only results concerning microquasars which may be observed by a 
telescope located in the Capo Passero site (i.e. with declination 5^ + 54°). For 
the most interesting steady microquasars the sensitivity is ~ 5- 10~ n erg/cm 2 
s, see Fig. 3. 
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Fig. 2. NEMO-km 3 sensitivity to neutrinos from microquasars versus source dec- 
lination, for a livetime of 1 year. The worsening of the sensitivity with increasing 
declination is due to the decrease of the source visibility. 



In Tab. 2 are summarized the detector sensitivities for microquasar LS 5039, 
according to the model of Aharonian et al. [8]. The simulations refer to a 
power-law neutrino spectrum, dn u /de u oc e~ r , with energy cutoff e™ ax = 10 
TeV and 100 TeV, and T = 1.5 and 2, respectively. For the four combinations 
of parameters T and e™ ax , an average neutrino energy flux fl h ~ 10~ 10 erg/cm 2 
s (e v > 0.1 TeV) was simulated. 
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Fig. 3. NEMO expected sensitivity for microquasar SS433 for the applied $™ ax and 
rbi n values as a function of the quality cut A C ut- 

6 Expected number of microquasar events 

In Tab. 3 are given the number of selected neutrino events from each micro- 
quasar, applying the cuts reported in Tab. 1 and according to the neutrino 
fluxes given by Distefano et al. [6]. The results refer to an integration time 
At equal to the duration of the considered burst for the transient sources and 
to one year for the steady sources. In the same table is given the background 
(atmospheric neutrinos + muons) in 1 year of data taking. In this analysis, 
it is assumed that transient sources cause one burst per year, i.e. the number 
of source events produced in the interval At is relative to 1 year observation 
time. 

In order to estimate the event rates for non-persistent sources it is crucial 
to know their duty cycle. Some of these sources have a periodic bursting 
activity: Circinus X-l has a period of 16.59 days [30]. This means therefore 
that we expect about 1.5 events per year. Other transient sources show a 
stochastic bursting activity. For such cases it is difficult to give an estimate of 
the expected event rate. For example, during 1994 GRO J1655-40 had three 
radio flares, each lasting 6 days [31]; during the same year GRS 1915+105 
emitted 4 bursts [32]. A recent study [33] has shown that microquasars GRS 
1915+105, Cygnus X-3 and Scorpius X-l are in flaring mode 21, 10 and 3 
percent of the time, respectively. The possibility to integrate over more then 
one burst could therefore help to detect neutrinos from microquasars. 

The search for neutrino events in coincidence with microquasar radio outbursts 
could be a tool to reject atmospheric background, restricting the analysis pe- 
riod to the flare duration At. Such an analysis technique, already used by 
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Table 1 

Detector sensitivity to neutrinos from micro quasars: the sensitivity /^go is 
calculated for an e~ 2 neutrino spectrum in the energy range 1 - 100 TeV, for a 
detector livetime of 1 year. The corresponding values of $™ ax , A cut and r& iri and the 
source declination 8 are also given. 



Source name 


r ( de §) 


A-cut 


rbin (deg) 


^,90 (erg/cm 2 s) 


5 (deg) 


Steady Sources 


LS 5039 


101 


-7.3 


0.9 


6.5 


10" 


ii 


-14.85 


Scorpius X-l 


104 


-7.7 


0.7 


5.8 


10" 


-11 


-15.64 


SS433 


115 


-8.0 


0.8 


5.7 


10" 


-11 


4.98 


GX 339-4 


96 


-7.4 


0.5 


4.7 


10" 


-11 


-48.79 


Cygnus X-l 


103 


-7.5 


0.7 


9.0 


10" 


-11 


35.20 


Bursting Sources 


XTE J 1748-288 


102 


-7.6 


0.9 


5.4 


10" 


-11 


-28.47 


Cygnus X-3 


101 


-7.3 


0.8 


1.1 


10" 


-10 


40.95 


GRO J1655-40 


101 


-7.4 


0.7 


5.2 


10" 


-11 


-39.85 


GRS 1915+105 


100 


-7.4 


0.8 


7.4 


10" 


-11 


10.86 


Circinus X-l 


90 


-7.3 


0.9 


4.2 


10" 


-11 


-56.99 


XTE J1550-564 


90 


-7.1 


0.9 


4.4 


10" 


-11 


-56.48 


V4641 Sgr 


102 


-7.4 


0.9 


5.6 


10" 


-11 


-25.43 


GS 1354-64 


90 


-7.5 


1.0 


3.8 


10" 


-11 


-64.73 


GRO J0422+32 


103 


-7.5 


0.8 


8.7 


10' 


-11 


32.91 


XTE J1118+480 


102 


-7.5 


0.7 


1.1 


lO- 


-10 


48.05 



Table 2 

Detector sensitivity flux to neutrinos from the steady source LS 5039: 

the sensitivity f u ,9o is calculated for a e~ r neutrino spectrum in the energy range 
from 0.1 TeV up to e™ ax and for a detector live time of 1 year. The corresponding 
values of •Q™ ax , A cut and run are also given. 



r 


e™ ax (TeV) 


tf™ ax (deg) 


Acut 


rbin (deg) 


fz/,90 (erg/cm 2 s) 


1.5 


10 


101 


-7.3 


1.0 


1.5 • 10~ 10 


1.5 


100 


101 


-7.3 


0.9 


5.2 • lO" 11 


2.0 


10 


97 


-7.4 


1.0 


2.7- 10- 10 


2.0 


100 


101 


-7.3 


0.9 


9.6 • 10~ n 



AMANDA [34], can improve the detector sensitivity to neutrinos from tran- 
sient sources. Referring to the bursts considered in Tab. 3 and integrating over 
the time interval At of the bursts, we expect an average background of about 
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10~ 3 events (muons) per burst. Summing on all the bursting sources in Tab. 
3, we count ~ 0.04 background events, which requires about 5 source events 
for a 5(7 level detection with a 70% probability [35]. Tab. 3 shows that we 
expect 3.4 + 9.0 events in the case of a burst from each of the bursting micro- 
quasars. Therefore, a cumulative analysis could provide a possible detection 
of microquasar neutrinos. 

For the microquasar LS 5039 we considered the flux predicted by Aharonian 
et al. [8]. The expected number of selected events is given in Tab. 4; the 
comparison with the atmospheric background shows that an evidence could 
be expected in a few years of data taking. 

In Tab. 5 we quote the expected number of events from microquasar Cygnus 
X-3, according to the Bednarek model [7]. Since the model predicts a neutrino 
flux with a spectral index close to 2 and since the main event contribution is 
in the energy range 1 - 100 TeV, we used the same event selection parameters 
quoted in Tab. 1. In the case of k — 2.5 and ^™ ax = 10 7 the signal is lower 
than the background. In the other cases, the number of events exceeds the 
expected background. The signal could be therefore detected in a few years of 
data taking, especially in the case of k = 2 and YA ax — 10 7 with 2.2 events 
per year. 



7 Conclusions 

The possibility to detect TeV neutrinos from Galactic microquasars with the 
proposed NEMO-km 3 underwater Cerenkov neutrino telescope has been in- 
vestigated. 

A Monte Carlo was carried out to simulate the expected neutrino-induced 
muon fluxes produced by microquasars and by atmospheric neutrinos. The ex- 
pected atmospheric muon background was also simulated. Muon tracks were 
propagated inside the detector using the ANTARES simulation tools, taking 
into account the water optical parameters measured in the Capo Passero site. 
Muon tracks were then reconstructed and the number of neutrino induced 
muons expected from microquasars was determined, according to different 
theoretical neutrino flux predictions quoted in the literature. These results 
were compared with the number of expected background events from atmo- 
spheric neutrinos and muons. We applied an event selection in order to reject 
the atmospheric background. The applied event selection is a combination of 
different criteria, chosen optimizing the detector sensitivity to neutrino fluxes 
coming from each microquasar. Eventually we calculated the expected number 
of events surviving the selection referring to 1 year of data taking. 
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Table 3 

Expected number of neutrino induced muons from the Levinson and 
Waxman microquasar model: N™ is the number of selected muons from each 
microquasar expected from the theoretical neutrino energy flux fl h quoted by Dis- 
tefano et al. [6], during the time interval At. We also report the expected number 
of atmospheric background events b surviving the event selection and expected in 1 
year of data taking. 



Source name 


At (days) 


fth ( er g/ cm 2 


s) 


NT 


b 


Steady Sources 


LS 5039 


365 


1.69-10- 12 




0.1 


0.1 


Scorpius X-l 


365 


6.48-10~ 12 




0.2 


0.1 


SS433 


365 


1.72-10~ 9 




76.0 


0.1 


GX 339-4 


365 


1.26-10- 9 




68.0 


0.1 


Cygnus X-l 


365 


1.88-10- 11 




0.5 


0.1 


Bursting Sources 


XTE J1748-288 


20 


3.07-10- 10 




0.8 


0.3 


Cygnus X-3 


3 


4.02-10~ 9 




0.8 


0.1 


GRO J1655-40 


6 


7.37-10- 10 




0.6 


0.1 


GRS 1915+105 


6 


2.10-10" 10 




0.1 


< 0.1 


Circinus X-l 


4 


1.22-10- 10 




0.1 


0.1 


XTE J1550-564 


5 


2.00-10- 11 




< 0.1 


< 0.1 


V4641 Sgr 


0.3 


2.25-10~ 10 ^-3.25 


• io- 8 


< O.l-rl.4 


0.1 


GS 1354-64 


2.8 


1.88-10- 11 




< 0.1 


0.1 


GRO J0422+32 


l-i-20 


2.5M0- 10 




< O.l-rO.4 


0.1 


XTE J1118+480 


30^150 


5.02-10- 10 




1.0-5-4.8 


0.2 



We made a comparison of our sensitivity to microquasar/point-like sources 
with the analysis made by IceCube [35,36]. For a 3 years exposure the sensi- 
tivities for a e~ 2 energy dependence are ~ 10~ 9 and ~ 2 • 10~ 9 GeV/cm 2 s for 
NEMO and IceCube, respectively. 



Our results show that, assuming reasonable scenarios for TeV neutrino pro- 
duction, the proposed NEMO telescope could identify microquasars in a few 
years of data taking, with the possibility of a 5a level detection; otherwise 
it would strongly constrain the neutrino production models and the source 
parameters. 
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Table 4 

Expected number of neutrino induced muons from microquasar LS 5039 
(Aharonian et al. model [8]): we give the number of selected events N™ from 
LS 5039 compared to the atmospheric background events b surviving the event 
selection, during an interval time of 1 year. 



LS 5039 



r 


e™ ax (TeV) 


N m 


b 


1.5 


10 


1.7 


0.2 


1.5 


100 


4.9 


0.1 


2.0 


10 


1.0 


0.3 


2.0 


100 


2.6 


0.1 



Table 5 

Expected number of neutrino induced events from microquasar Cygnus 
X-3 (Bednarek model [7]): N™ is the number of selected muons from Cygnus 
X-3 during an interval time of 1 year. We also report the number of atmospheric 
background muons b surviving the event selection during the same interval time. 



Cygnus X-3 


REG. II 


REG. Ill 






N m 

A 4 


N m 

M 


b 


2.0, 10 6 


0.4 


0.4 


0.1 


2.0, 10 7 


1.1 


1.1 


0.1 


2.5, 10 7 


< 0.1 


0.1 


0.1 


mono-energetic 


< 0.1 


0.7 


0.1 
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